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Nonequilibrium processes in semiconductors are considered with highly nonuniform initial densi- 
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ties of charge carriers. It is shown that there exist such distributions of charge densities under which 
the electric current through a sample displays quite abnormal behaviour flowing against the applied 
voltage. The appearance of this negative electric current is a transient phenomenon occurring at the 
initial stage of the process. After this anomalous negative fluctuation, the electric current becomes 
normal, i.e. positive as soon as the charge density becomes more uniform. Several possibilities for 
the practical usage of this effect are suggested. 
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^j- , I. INTRODUCTION 

o , 

The study of electric processes in semiconductor materials plays an important role in understanding the physics 
of semiconductor devices as well as in their design and development [1,2]. One of the most difficult problems is the 
consideration of strongly nonequilibrium effects in essentially nonuniform semiconductors. At the same time the latter 
can display quite interesting specific features caused by nonuniform distributions of charge carriers [3-6]. 

For instance, electric current through a semiconductor device can display rather abnormal behaviour, with transient 
""^J ■ fluctuations corresponding to the flow of the current against the applied voltage [5,6]. In the short communications 
[5,6] a model case of a unipolar semiconductor was considered, with numerical analysis not including the relaxation 
parameters and diffusion coefficients. 

The aim of the present paper is to study nonequilibrium processes, with a strongly nonuniform initial distribution 
of charge carriers, for realistic semiconductor materials. We consider the general case of a semiconductor with two 
kinds of charge carriers, positive and negative. Numerical analysis takes account of relaxation and diffusion effects. 
The total current through a semiconductor device is considered, together with the currents across the left and right 
surfaces of this device. And also the role of the generation-recombination noise is analysed. 

In Sec. 2 we collect the main equations defining the problem and needed for the following analysis. An approximate 
analytical solution of these equations is presented in Sec. 3, which permits to show explicitly the motion of charge 
carriers under an applied voltage. The conditions for the occurrence of the negative electric current, directed against 
this applied voltage, are derived in Sec. 4. This anomalous transient effect is illustrated by numerical solutions. Sec. 5 
contains conclusion and discussion with several suggestions for the possible practical usage of the considered effect. 
Appendix contains the proof that the approximate regular solution asymptotically coincides with the regular stable 
solution of the exact equations. 



II. BASIC EQUATIONS 

The charge carriers are characterized by the densities p\ > and p2 < that are functions of the space 
position r and time t , i.e. pi = pi(r ,<) with i = 1,2 . Carrier transport can be described in terms of a 
semiclassical approach called the drift -diffusion approximation which is the basis of the majority of semiconductor 
device models [1,2]. This approach is, of course, phenomenological since it does not concern the microscopic derivation 
of the used parameters, such as mobilities or relaxation widths. The values of these parameters depend on a number 
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of different underlying causes. For instance, the lattice structure of the considered semiconductor is important for 
defining the values of these parameters. Thus, the scattering of carriers on phonons influences both mobilities as well 
as damping parameters. However, the calculation of such parameters is a separate problem that is not the aim of 
the present paper. In the semiclassical approach all parameters are assumed to be given a priori. With the given 
phenomenological parameters, the drift-diffusion approximation is known to give a very good description of realistic 
semiconductor devices, which explains why this approximation is so widely used [1,2]. The first set of equations in 
this approach consists of the set of continuity equations 

at n 

for each kind of carriers. The relaxation term, with the relaxation time , is taken in the simplest form since, as will 
be clear from what follows, it does not play an essential role in transient processes occurring at times t <C n . The 
right-hand side of (1) is the generation-recombination noise [7] which is always present in semiconductor devices. Not 
to overcomplicate the problem we do not include other types of noise [7-9] assuming that they are of less importance. 
Another equation is a Maxwell equation 



with the dielectric permittivity e . 
The electric-current density in (1) is 



eV ■ E= MPi + Pi) , (2) 



3 i= P-iPi E —Di V Pi, (3) 



where ji\ > and /j,2 < are the carrier mobilities and Di = (/li/e^ksT are the diffusion coefficients in which 
ei > and e 2 < are the carrier charges, fcs is the Boltzmann constant, and T is temperature. The first and 
second terms in (3) are the drift and diffusion current densities, respectively. Adding here the displacement current, 
one has the total current density 

-? t -? e dE 

3tot=3i + J 2 +^^f ■ (4) 
If the device is biased with an externally applied voltage Vq along a path I , then 

E(r,t)d7=V Q . (5) 



For concreteness, we consider a positive bias, that is, Vb > . It is assumed that metal contacts supplying the 
external voltage are nondamaging, that is, do not induce in their vicinity incubation effects. 

Consider a plane device of the width L and area A . Then instead of r we are to deal with one space variable 

x e [o, L] . 

An important characteristic is the transit time 

where \i = min{/ii, I/12I} ■ Usually, fi\ < l^l • This time is to be compared with the relaxation times Tj . It is often 
more convenient to deal with inverse times called widths. For instance, 

7« = -, (7) 

n 

with i = 1, 2 , define the relaxation widths. 

Another convenience is to deal with dimcnsionlcss quantities, which will be done in what follows. To return to 
dimensional quantities, we shall imply that x is measured in units of L ; t and Tj , in units of r ; and other 
physical quantities, in the corresponding units listed below: 

D a = M V , E a = Q = eAE , 
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PO=-TT, SO = — , Jo = • (O) 

AL T r 



For the case considered, Eq.(l) reduces to 



and Eq.(2) becomes 

dE 

— =4n( Pl + P2 ), (10) 

where 

< x < 1, t>0. (11) 

Condition (5) for an applied voltage reads 



f E(x,t)d. 
Jo 



\x = 1. (12) 

These equations are to be supplemented by initial conditions 

Pi (x,0) = fi(x) (i = 1,2). (13) 
Electric field can be expressed, from Eq.(10), as a functional 



Q(x, t) — I Q(x, t)da 
Jo 



(14) 



E(x,t) = 1 + 4tt 
of the charge densities, so that 

Q(x,t)= f [p 1 (x' 1 t) + p 2 (x' 1 t)}dx'. (15) 
Jo 

And the total density of current (4) can be written as 

3 „ = tufl E-n^ +mP2 E-D 2 ^ + ±^, (I.) 

where j to t = jtot(x,t) . 

The quantities that can be measured and that we are going to study in what follows are the current across the left 
surface 

J(0,t) = j tot (0,t) , (17) 

the current across the right surface 

J(l,t) = j tot (l,t) , (18) 

and the total current through the device 

J(t) = f jtot(x,t)dx . (19) 
Jo 

The latter, using Eqs. (16) and (12), can be presented as 

J{t) = \ [pipi(x,t) + H2p 2 (x,t)}E(x,t)dx+ 
Jo 

+ D 1 [p 1 (0,t)-p 1 (l,t)]+D 2 [p 2 (0,t)-p 2 (l,t)} . (20) 

Our aim is to study the peculiarities in the time dependence of the electric current when the initial conditions (13) 
correspond to a strongly nonuniform charge distribution. Such nonuniform distributions can be prepared in different 
ways. For example, one can organize a nonuniform distribution in the process of growing of a semiconductor sample. 
Another way is to irradiate semiconductor by narrow laser beams [10]. One more possibility is by forming heavily 
doped layers by ion irradiation [11]. This method makes it possible to form narrow layers of positive carriers with a 
density of 10 20 cm~ 3 . 



3 



III. CARRIER DENSITIES 



To understand better the physics of processes resulting from a nonuniform initial distribution of charge carriers, 
it would be useful to find an analytical, though approximate, solution to the system of equations (9) and (10). This 
can be done by means of the method of scale separation [12,13], whose mathematical foundation is based on the 
Krylov- Bogolubov averaging method [14]. 

The first step in the method of scale separation [12,13] is to classify the solutions onto fast and slow. In our case 
this can be done as follows. The electric field, as is seen from Eq. (14), is the functional of the charge densities, 
averaging the latter over the space variable x . This results in that E varies in space slower than pi . On the other 
hand, the voltage integral (12) shows that the electric field, being averaged over space, does not depend on time. This 
means that E can be treated as a slow function in time. Therefore, the electric field E can be regarded as a slow 
solution, as compared to the charge densities pi , with respect to both space and time. This permits us to consider 
the equation (9) for a fast solution pi keeping there E as a space-time quasi-integral. After solving Eq. (9), the 
found pi is to be substituted into Eq. (14) giving an equation for E which can be solved iteratively. 

In solving Eq. (9), it is convenient to continue pi outside the region of x G (0, 1) by defining pi as zero for 
x < and x > 1 . Then we may invoke the Fourier transforms with respect to x . Finally, we obtain 

Pi {x,t)=pY a {x,t)+pT n {x,t); (21) 

the first term being the regular solution 

/+oo 
Giix-x'^Mx^dx' (22) 
-oo 

induced by the initial condition (13), while the second term being the random solution 

rt p + oo 

p r l an (x,t)= / Gi(x-x',t-t%(x',t')dx'dt' (23) 

JO J-oo 

generated by the noise. The Green function in Eqs. (22) and (23) is 

1 [ + °° 

Gi(x,t) = — / exp{ikx ~ ibJi(k)t}dk (24) 

with the spectrum 

uoi{k) = p,iEk — iDik 2 — iji. (25) 



Function (24) has the properties 



r+oo 

Gi(x,0) = S(x), I G i (x,t)dx^e-^ t . 
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In the case considered, the integration in (24) can be realized explicitly resulting in 

^•^^M-^r 1 "4 (26) 

As the initial condition in Eq. (13) it is reasonable to accept the physically realistic case of the Gaussian distribution 

fi(x) = — ex P ^ — r ' ( 27 ) 



in which < dj < 1 and 

Qi = J fi(x)dx, Z, = J cxp ^ jeb- 

With the initial condition (27), the regular solution (22) becomes 
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p^(x,t) = ' cxp^ v ,7 l; -Jit}. (28) 

The regular and random solutions satisfy the initial conditions 

p r i eg (x,0) = f i (x), p™ n (x,0) = 0. (29) 
The regular solution, as time increases, moves with the velocity piE , becomes wider and smaller, so that 

lim p r C g {x,t) = 0. (30) 

t — >oo 

The behaviour of the random solution depends on that of the noise. It is customary to treat the latter as the white 
noise with the averaging properties 

(Si(x,t))=0, ^ i (x,t)^(x',t , )}=- fij S(x-x')S(t-t'). (31) 

Accepting (31), one has 

(p™ n (x,t))=0, (32) 

and, consequently, 

lim (pi(x,t)) = 0. (33) 

t— >oo 

Since the electric field, according to (14), is a linear functional in pi , we find 

lim {E(x,t)) = 1. (34) 

t — >oo 

Although the limiting values (33) and (34) have been obtained by analysing the approximate solutions, it is possible 
to show (see Appendix) that the limits (33) and (34) are stable stationary solutions of the exact equations. 
The random solution (23) influences the electric current (2) through the correlator 

(pl an (x,t)p r j an (x',t))= lij f G t3 {x - x 1 ,t')dt' , (35) 

Jo 

in which 

/+oo 
G l {x - x",t)Gj{x' - x",t)dx". (36) 
-oo 

With the Green function (26), this gives 

^ , x 1 f (x - PiEt + a.Etf N 1 

G " M = v?IB=rsr x,> {- m+w - H + ^1 ■ (37) 

At large time, Eq.(37) decays by the law 

Gij (x, t) ~ 1 = cxp (-7 e //t) , (38) 

Z^TT{L>i + L)j)t 



as t — > oo , with the effective attenuation 



(p*-Pj) 2 E 2 
leff ~ i(Di + Dj) + 7<+7 ^ 



Consequently, the correlator (35) tends to a time constant. 
At small time, one has 



G ^^ 2 MD ] +D y A-Mo^Dj t ^ (39 » 



as t — ► . Therefore, the correlator (35) behaves as 



( P r(x,tWr(x',t)) - ^£== ex p(- -^ ^ ^ - (40, 



2 v /tt(A + D/) 4(A + !>,•)«;' 

when i — > . Equation (40) shows that the influence of noise at small times is exponentially suppressed. This 
conclusion is of high importance for the following analysis. 
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IV. ELECTRIC CURRENT 



We have now enough information about the behaviour of the system in order to answer the question: Is it possible 
that a negative electric current could appear, directed against the applied voltage? 

The first evident necessary condition for such a possibility is the space nonuniformity of the charge densities. Really 
if pi(x,t) is uniform in x , then from Eq. (20) it follows immediately that the electric current is the positively 
defined quantity p\p\ + P2P2 > . 

The properties of the carrier densities, studied in the previous section, are such that, even, if at the initial time 
Pi(x,t) is nonuniform in space, it tends to become uniform with time. Consequently, if a negative electric current 
would appear, this could happen only at the initial stage of the process, when tel. 

In this way, the occurrence of negative electric current, if any, can arise only as a principally transient effect, 
when the charge densities are yet nonuniform. The processes of diffusion and relaxation need some time to make 
Pi uniform. Therefore, there always can be found such a time f«l when diffusion and relaxation are yet not 
important. But these processes shorten the time of a negative-current fluctuation, if it appears. Thus, the conditions 
favoring the longer lifetime of such a fluctuation would be and 7^ <C 1 . 

The influence of noise, according to Eq. (40), is exponentially small at the initial stage. Thus, even a strong noise 
would not kill the effect, although it, of course, would shorten the negative-fluctuation lifetime. So, the condition 
favoring the longer lifetime is weak noise, when 7^ -C 1 . 

After understanding the necessary and favoring conditions for the transient effect of a negative-current fluctuation, 
let us elucidate sufficient conditions for the inequality 

J(t) < (41) 

occurring at t <C 1 . As a limiting case we may take t = and the maximally nonuniform initial density 

Pi (x,Q) = U{x) = Qi&{x - <n) (42) 

following from condition (27) under bi — > . Then from expression (20) we readily get 

J(0) - fJLiQiE{ai,0) + p 2 Q2E(a 2 ,0). (43) 

This emphasizes once again that for such a nonuniform initial charge density the diffusion, relaxation, and noise do 
not influence much the value of the electric current J(0) . The corresponding electric field, defined by Eq. (14), is 

E(x, 0) = 1 + 47rQi[oi - 0(oi - x)] + 4irQ 2 [a 2 - 9(a 2 - x)], (44) 

where Q(x) is the unit step function. 

Combining Eqs. (41), (43), and (44), we have the inequality 



+ 



P1Q1 |l + 4ttQi ^ai - +4ttQ 2 [a 2 - 9(a 2 - ai)] j 

+ P2Q2 |l + 4ttQ 2 [a 2 - +4ttQi [ch - 0(ai - a 2 )]| < (45) 

for the parameters of the system allowing the appearance of a negative current. 

There can be a number of different cases satisfying inequality (45). To show that such situations do really exist, 
consider a particular example when a\ = a 2 = a . Then Eq.(41) reduces to 

(piQi+p 2 Q2)E(a,0) <0. (46) 

From equality (44) we get 



£(a,0) = 1 + 4ttQ (a- {Q = Qi + Q 2 )- 

Recall that pi and Qi are positive, while p 2 and Q2 are negative; so that piQi > . Thence, inequality (46) 
can be hold only if 

£7(o,0)<0. (47) 
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As follows from solution (28), the quantity y^E plays the role of the effective velocity of motion for the corresponding 
charge packet. In the case of inequality (47), we have [i\E < and > . This means that the positive carriers 
effectively move against the applied voltage; and the negative carriers, along the latter; that is, they move oppositely 
to what one would expect. Hence, the negative electric current is related to the anomalous drift of charge carriers. 
Substituting into Eq. (47) the electric field, wc find 

4ttQ ( I - a) > 1. (48) 



2 

Depending on whether Q is positive or negative, inequality (48) yields 



a> l + ^Q\ (Q<0) ' (49) 
Taking also into account that < a < 1 , wc obtain from inequalities (49) the condition 

\Q\ > (50) 

Equations (49) and (50) are sufficient conditions for the appearance of a negative electric current at the initial stage 
of the process. Similarly, it is easy to show that the currents (17) and (18) can also become negative in the transient 
regime. 

To study in more detail the behaviour of the electric current as a function of time, we have solved Eqs.(9) and (10) 
numerically. In accordance with the above analysis, the case favoring the negative-current fluctuation is accepted, 
when 7y <C 1 . The initial conditions are given by the Gaussian form (27). The voltage integral (12) plays the 
role of the boundary condition for the electric field. For the charge densities one may take the Neumann or Dirichlct 
boundary conditions [1]. We have tried both and found that the general picture does not change much, with the only 
difference that the calculational procedure is less stable for the Dirichlct conditions. To achieve the best stability, we 
opted for the Neumann boundary conditions. 

For the characteristic parameters we accept the values typical of semiconductors [1,2], such as Si . Then the 
diffusion coefficients are D\ ~ 10 cm 2 /s, D 2 ~ 30 cm 2 j s. The mobility of positive carriers \x\ ~ 500 cm 2 /Vs for the 
average concentration 10 13 — 10 14 cto~ 3 and \x x ~ 200 cm 2 /Vs for the concentration 10 17 — 10 18 cm~ 3 . The mobility 
of electrons /j,2 ~ 1500 for the average concentration 10 13 — 10 14 and y,2 ~ 400 for the concentration 10 16 — 10 18 . 
The recombination time t\ <~ t 2 ~ 10~ 12 — 10~ 10 s , hence the relaxation width 71 <~ 72 ~ 10 10 — 10 12 s _1 . We 
consider a plane device of the size A <~ 1 cm 2 , L ~ 0.1 — 1 cm , with the applied voltage Vo ~ 10 3 — 10 5 V . 
For the calibration parameters in Eq.(8) we get fj, ~ 10 3 cm 2 /Vs, V ~ 10 3 - 10 5 V, D ~ 10 6 - 10 8 cm 2 /s, E ~ 
10 3 - 10 5 Vycm, Q ~ 10 3 - 10 6 V cm, p ~ 10 3 - W 7 V/cm 2 , J - 10 8 - 10 15 y cm 2 /s. The transit time (6) is 
T - 10~ 6 - 10~ 9 s. 

The results of our numerical calculations are presented in the Figs.l to 6, where we show the time dependence of the 
electric current across the boundaries as well as the behaviour of the total current through the device. These figures 
demonstrate that the appearance of negative electric current is really a transient effect occurring at dimcnsionlcss 
times id, which in dimensional units means that t <C tq . All values in the figure captions are given in dimensionlcss 
units employing the calibration parameters from Eq. (8). Also, for shortness, we write a\= a and b\ = b. 

Fig.l shows that the electric current through the left boundary of the semiconductor sample, through its right 
boundary, and the total electric current are different. This difference is not merely quantative but can be qualitative, 
so that the negative current may happen at the right surface and on average through the sample, but may be absent 
on the left surface. Such a difference depends on semiconductor characteristics as mobilities and relaxation widths. 
This suggests the possibility of employing the principal difference in the behaviour of the currents for extracting 
information on the semiconductor characteristics. For example, in Fig. 2 it is seen that changing the electron mobility 
mainly influences the current across the left surface. The transient negative current becomes more pronounced when 
increasing the absolute value of the electron mobility, as is seen in Fig. 3. The lifetime of the negative-current 
fluctuation strongly depends on the relaxation width, which is illustrated in Fig. 4. Increasing the relaxation width 
shortens the fluctuation lifetime. Figure 5 demonstrates the role of the total initial charge on the occurrence of the 
negative current, and Fig. 6 shows the role of the initial distribution of charge carriers. The importance of special 
conditions for the initial charge and its location has been discussed in detail above. The amplitude of the transient 
negative-current fluctuation becomes smaller when the charge layer at initial time is shifted farther from the left 
surface of the semiconductor sample. 
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V. DISCUSSION 



We have considered electric processes in nonequilibrium nonuniform semiconductors. The transport equations are 
taken in the standard drift diffusion approximation that is widely used for describing realistic semiconductor devices. 
Both analytical and numerical solutions of these equations are accomplished. 

It is shown that under special circumstances an unusual transient phenomenon appears displaying negative electric 
current. The necessary condition for such an anomalous current is nonuniformity of carrier densities at the initial 
stage. A general sufficient condition (45) is derived and its particular forms (49) and (50) are analysed in detail. We 
studied the influence of diffusion, relaxation, and of generation-recombination noise and showed that these processes, 
even being strong, do not destroy the effect although may shorten the lifetime of a negative-current fluctuation. 
Therefore such an anomalous electric current can really be observed in semiconductor devices. 

An important physical question is how one could use the considered effect for practical applications. Several 
possibilities of using this effect can be suggested: 

(i) The appearance of the transient negative-current fluctuation is rather sensitive to the characteristic parameters 
of semiconductor, such as the carrier mobilities and relaxation widths. Therefore one could use the observation of the 
current fluctuation for defining these parameters. This could be done in the following way. Assume that we know all 
parameters except one, say a mobility or a relaxation coefficient. Comparing the time-dependence of the measured 
current with that of the calculated one, we may try to find such a value of the sought parameter that the measured 
and the calculated behaviour of the electric currents be as close as possible to each other, in the optimal case, be 
almost coinciding. Then the so fitted quantity would give the value of the sought parameter. 

(ii) When the layer of charge carriers is formed by irradiating semiconductor with an ion beam, the stuck ions are 
distributed approximately in the Gaussian law centered at the mean free path of the ions. A necessary condition for the 
occurrence of the negative current is that the charge layer is located at a particular distance from the semiconductor 
surface, as e.g. in Eq. (49). Thence, this effect is very sensitive to the initial location of charge carriers and, thus, 
could be used for measuring the mean free path of ions in specific semiconductor materials. 

(iii) The value of the total charge in the initial nonuniform layer is also crucially important for the occurrence of the 
negative current, as is seen from Eqs. (48)— (50). Hence, studying this current, we could measure the initial charge. 
The latter may be unknown when the initial distribution of charge carriers is formed by irradiating semiconductor 
with narrow laser beams whose influence on the generation of carriers in not precisely known. 

(iv) Semiconductor devices often work in the close vicinity of radiation sources, such as atomic reactors, or under 
the influence of other strong radiation, as cosmic rays. In the presence of radiation, the functioning of semiconductor 
devices can be drastically disturbed because of the arising carrier nonuniformities. This can lead not only to the 
malfunctioning of semiconductor devices but even to dramatic accidents. In order to prevent from these, one could 
employ controlling schemes reacting to the appearance of the negative current, signaling by this that the level of the 
carrier nonuniformity induced by irradiation has become dangerous. 

(v) As we have shown, the generation-recombination noise does not destroy the effect of the negative-current 
fluctuation. However, there exist other types of noise [7-9] whose influence on the supression of this effect can be 
different. Therefore, analysing the peculiarity of the electric negative-current fluctuation, one could judge what kind 
of noise dominates the process in the studied semiconductor. 

It is certainly not possible to enumerate all feasible applications of the considered effect. But we hope that the 
examples listed above do demonstrate that the specific unusual features of the negative-current transient effect could 
provide us several interesting physical applications. Three types of such applications are, generally admissible. One 
type is for investigating the characteristics of semiconductor materials. Another type can be used for studying the 
properties of irradiating beams. And, finally, this effect can be employed for the practical purpose of creating special 
controlling instruments. 
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Appendix. Stationary Solution 



Here we prove that Eqs.(33) and (34) are the stable stationary solutions of the exact equations (9) and (10) under 
the voltage condition (12). To this end, because of equality (32), it is sufficient to prove that 

lim p\ e \x, t) = 0, lim E reg (x, t) = 1, 

t — >oo t — »oo 

where E reg implies the functional (14) including the dependence on only p^ 9 . In what follows we shall write for 
brevity pi and E keeping in mind p\ e9 and E re9 . The proof will be based on the method of multipliers [15]. 
Write Eq.(9) in the form 

d 

—Pi(x,t) = Vi(x,p,t), 



with the velocity field 



Define the multiplier matrix 



and the Jacobian matrix 



d 2 pi 



Mij(x,x',t) 



Lij(x,x',p,t) = 



Spi(x,t) 
6 Pj (x',0) 

5vi(x,p,t) 



Spj(x',t) 

The latter, with the given velocity field, consists of the elements 



f d 2 d \ 

L u (x,x ,p,t) = ( D t — - p t E— - 7i S(x-x')- 



-^pi(2pi + pj)S(x - x') - p. 



dpi SE(x, t) 
dx 5pi(x', t) 



Lij(x, x', p, t) = -4:TTfj,iPiS(x - x') - p, 



dp t SE(x,t) 



dx 5pj(x', t) ' 

where i ^ j . Here, the variational derivative of the electric field can be found from expression (14) giving 



6E(x,t) 
5pi(x',t) 



47T 



6Q(x,t) 
Spi(x',t) 



+ x' - 1 



5Q(x,t) 
Spi(x',t) 



e(x-x'). 



Varying the evolution equation for pi , we obtain the equation 



^-M ij {x,x',t) = y2 I L ik (x,x",p,t)M kj (x",x',t)dx" 
at i. Jo 



for the multiplier matrix, with the initial condition 

Mij(x, x' , 0) = Sijd(x — x'), 

following from the variation of condition (13). 

For the case pi — and E = 1 , the Jacobian matrix is 



/ d 2 d \ 
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This leads to the equation 



dM t1 ( d 2 d 



dt \ dx 2 dx 
for the multiplier matrix. From the latter equation, invoking the Fourier transform 



Muixrfj) = — / M l3 {k,ty k{x ~ x Uk, 

27T J-oo 



we find 

Mij(k,t) = 6 ij exp{-(ifi i k + D, L k 2 +^i)t] 

As is seen, 

\Mij(k,t)\ < 1 

for all k e (—00, +00) and t > . Hence the motion in the vicinity of the stationary solutions pi = and E = 1 
is stable. It is also asymptotically stable since 



t— >oo 

for all k e (— 00, +00) . This completes the proof. 



lim Mij(k,t) = 
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Figure Captions 



Fig.l The time behaviour of the electric current (17) across the left boundary (solid line), of the current (18) across 
the right boundary (long-dashed line), and of the total current (19) through the semiconductor device (short-dashed 
line). The used parameters, in dimcnsionless units, are Q\ = 1, Q 2 = —\, Mi = 1; M2 = — 10, 7i = 72 = 10, a — 

0. 25, b = 0.1, £>i = 1(T 3 , D 2 = 3D X . 

Fig. 2 The same as in Fig.l, but for [i 2 = —3 and a = 0.05. The smaller absolute value of the electron mobility 
influences mainly the current across the left boundary. 

Fig. 3 Electric current (19) through the semiconductor device as a function of time for different mobilities. The 
fixed parameters are Qi — 1, Q 2 — — 2 , 7i = 72 = 1, a = 0.25, b = 0.1, D\ = 10~ 3 , D 2 = 3D a , fii = 1. The curves 
correspond to the varying mobility: \i 2 = —10 (solid line), fi 2 = —5 (long-dashed line), and fi 2 = —3 (short-dashed 
line). Increasing the absolute value of the electron mobility makes the negative current more pronounced. 

Fig. 4 The same as in Fig. 3, except for [l 2 = —3 and for varying relaxation widths: 71 = 72 = 25 (solid line), 
71 = 72 = 10 (long-dashed line), and 71 = 72 = 1 (short-dashed line). Increasing relaxation width suppresses the 
negative current. 

Fig. 5 The time dependance of the electric current (19) for the following parameters: m = 1, fj, 2 = —3, 71 = 72 = 

1, a = 0.25, b = 0.1, Di = 10~ 3 , D 2 = 3Di, Q\ = 1, and for varying Q 2 = (solid line), Q 2 = -\ (long-dashed 
line), Q 2 = — \ (short-dashed line), Q 2 = — | (dotted line), and Q 2 = — 1 (dash-dotted line). 

Fig. 6 The same as in Fig. 5, but for the fixed Q 2 = — \ and for a varying location of the initial distribution of 
carriers: a = 0.1 (solid line), a = 0.2 (long-dashed line), a — 0.3 (short-dashed line), a = 0.4 (dotted line), and 
a = 0.5 (dash-dotted line). 
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